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ABSTRACT 


In  this  thesis,  the  amplitude  gain  of  received  quadrature  phase  shift  keying  (QPSK) 
signals  and  interferences  are  estimated  via  the  least  squares  error  method  in  order  to 
facilitate  the  implementation  of  reference-based  successive  interference  cancellation 
(RSIC).  In  the  scenario  considered  for  this  thesis,  a  system  of  multi-platform  receivers  is 
assumed  to  be  positioned  within  the  overlapping  coverage  areas  of  transmitting  base 
stations.  The  first  receiver  initiates  the  RSIC  technique  by  obtaining  a  demodulated 
reference  signal  and  forwarding  that  reference  to  the  second  receiver  that  separately 
collects  a  second  signal  corrupted  by  interference,  which  is  actually  a  scaled  version  of 
the  first  signal  with  an  unknown  amplitude  gain.  This  amplitude  gain  is  estimated  and 
applied  to  the  known  reference  signal  so  that  it  is  subtracted  from  the  collected  signal  at 
the  second  receiver.  The  process  continues  with  a  third  and  fourth  receiver  (or  potentially 
more  receivers)  until  the  final  desired  signal  is  separated  from  all  the  interferences. 
Finally,  the  accuracy  of  the  estimations  is  evaluated  and  the  symbol  error  rate 
performances  via  Monte  Carlo  simulations  for  QPSK  modulation  in  a  multi -platform 
system  are  presented. 
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EXECUTIVE  SUMMARY 


Wireless  communication  has  made  huge  advances  in  the  ability  to  transmit  and  receive  at 
high  data  rates  in  almost  every  populated  region  in  the  world.  The  amount  of  interference 
experienced  by  these  systems  may  seem  overwhelming  at  times,  but  innovations  in 
technology  and  modulation  continue  to  provide  the  world  with  reliable  exchange  of 
information.  Another  technique  exists  that  is  utilized  by  a  multi-platform  system  to 
remove  extraneous  interference  from  a  received  signal  so  that  only  the  desired  signal  of 
interest  (SOI)  remains. 

Successive  interference  cancellation  (SIC)  is  considered  a  member  of  the  family 
of  post-interference  cancellation  techniques  which  are  focused  on  the  cancellation  or 
suppression  of  interference  on  an  a-posteriori  basis  [1].  This  means  that  a  system  using 
the  low  complexity  SIC  method  must  first  collect  a  signal  in  order  to  then  “demodulate 
and/or  decode  desired  information,  and  use  this  information  along  with  channel  estimates 
to  cancel  received  interference  from  the  received  signal”  [2].  A  related  SIC  scheme 
known  as  reference-based  successive  interference  cancellation  (RSIC)  makes  use  of  a 
multi-platform  system  with  what  is  called  a  reference  signal  to  eliminate  interference  in 
received  signals. 

RSIC  makes  use  of  multiple  receivers  placed  in  favorable  locations  so  that  each 
can  collect  a  specific  combination  of  signals.  An  iterative  process  ensues  in  which  a 
receiver  in  one  area  successfully  collects  a  signal  (called  the  reference),  demodulates  and 
then  transmits  it  to  the  successive  receiver  in  the  next  area.  Reference  signals  are 
assumed  to  be  mostly  interference-free.  The  second  receiver  separately  collects  a 
different  signal  in  addition  to  interference,  which  is  a  scaled  version  of  the  first  signal 
(scaled  due  to  amplitude  gain  or  loss).  By  calculating  an  estimate  of  the  amplitude  gain 
and  applying  it  to  the  reference  signal  that  is  obtained  from  the  first  receiver,  the  second 
receiver  can  subtract  this  interference  from  the  overall  signal  it  has  received.  This 
effectively  provides  the  second  receiver  with  its  own  reference  signal  to  pass  on  (in 
addition  to  the  first  reference  signal)  to  the  next  receiver  in  the  system  and  thus 
continuing  the  RSIC  process. 
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Previous  work  related  to  RSIC  uses  this  described  technique  to  model  systems 
with  different  amplitude  gain  configurations  with  the  intention  of  observing  the 
performance  provided  by  the  interference  cancellation  [3].  For  that  study,  the  amplitude 
gain  parameter  values  are  chosen  arbitrarily  simply  to  examine  how  the  cancellation 
technique  functions  in  different  interference  environments.  The  purpose  of  this  thesis  is 
therefore  to  provide  a  method  to  estimate  the  values  of  the  gain  parameters  experienced 
by  a  multi-platform  system,  in  order  to  allow  the  RSIC  technique  to  successfully  remove 
extraneous  unwanted  interference  from  SOIs.  The  performance  and  accuracy  of  these 
estimates  will  also  be  analyzed. 

As  seen  in  Figure  1,  the  overlapping  base  station  (BS)  coverage  areas  contain 
sectors  where  different  sets  of  signals  coincide. 


Figure  1.  Configuration  of  multiple  BSs  with  overlapping  antenna  coverage  which  has 
different  signal  combinations  received  in  each  sector.  From  [3]. 

The  received  signals  from  each  sector  are  given  by 


y,  =5i+wp 

(1.1) 

y2  S2  ^'1,2^1  ^2  ’ 

(1.2) 

y3  —  s3  +  a2  3s0  +  al  3.s,  +  w3 , 

(1.3) 

^4  ^  ^3,4^3  ^  ^'2,4^2  ^  ^1,4^1  ^  ^4  ? 

(1.4) 
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where  am  n  is  the  amplitude  gain  (or  loss)  of  the  signal  from  BS  m  received  by  receiver  n, 

and  w  is  the  additive  white  Gaussian  noise  in  the  receiver.  The  transmitted  signals  s  are 
assumed  to  have  perfect  synchronization  and  the  same  energy  while  the  amplitude  gain  is 
assumed  to  be  real  and  positive.  Receivers  in  each  sector  are  assumed  to  only  have  the 
knowledge  of  the  signal  y  they  collect,  and  the  reference  signals  that  are  sent  to  them  by 
the  other  receivers  (i.e.,  receiver  two  (Rx2)  only  knows  y2  and  y  ;  Rx3  only  knows  y3 , 

s2  and  y ;  Rx4  only  knows  y4 ,  s3,  s2  and  y ).  If  for  example  s2  is  considered  the  SOI, 

then  in  order  to  conduct  RSIC,  the  receiver  placed  in  Sector  A  would  collect  (1.1), 
demodulate  the  information,  then  transmit  this  reference  signal  to  receiver  two.  At  Rx2, 
the  received  signal  y2  and  the  obtained  reference  signal  would  be  used  to  calculate  an 

estimate  for  the  amplitude  gain  ax  2  so  that  the  estimate  could  be  multiplied  to  the 
reference  and  then  subtracted  from  the  received  signal  y2 .  The  SOIs  from  each  sector  are 
therefore 
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(1.5) 

a  =  y2~ P\2sx  ~w2> 

(1.6) 

A  —  -A  —  02,3s 2  ~  0i,3Sl  —  'A  ’ 

(1.7) 

A  —  y*4  —  03,4S3  —  02,4s 2  ~  0\,4Sl  ~  W4’ 

(1.8) 

where  /?  is  the  amplitude  gain  estimate  for  the  respective  a  parameter  and  the  noise  w 
may  be  considered  negligible  for  high  signal-to-noise  ratios  (SNR). 

In  order  to  obtain  these  [3  parameters,  the  least  squares  error  (LSE)  method  from 
[4]  is  used  to  calculate  estimations  for  each  sector.  In  the  example  where  s2  is  the  SOI, 
Rxl  collects  y,  which  it  passes  through  a  maximum  likelihood  detector  (MLD)  [5]  to 
decide  on  what  information  is  received,  or 


Sj  =  dec(  y , )  =  dec(sx  +  w, ), 


(1.9) 


where  dec(  )  is  the  MLD,  y,  =  [y1[0]y1[l]...;y1[./V-l]f  and  s,.  =  [s,[0]s,[l]...s,.[Ar-l]f 
for  i  =  l,2.  Next,  the  decision  as  to  what  symbols  are  received  by  Rx2  is  made  based  on 
y2  alone,  such  that  we  have  an  initial  (pre )  calculation  for  s2,  or  simply  s2  where 

s2  =  dec(y2).  (1.10) 

This  pre  calculation  is  used  to  determine  an  estimate  for  ax  2 ,  which  is  defined  as  fdX2.  The 
value  for  f5x  2  is  now  used  to  determine  a  final  (post)  calculation  for  the  information  of  s, 
in  (1.6),  or 

s2=dec(y2-/3l2sl).  (1.11) 

If  the  concept  of  making  decisions  to  obtain  reference  signals  is  executed  on  a 
three-receiver  system,  by  utilizing  s3  =  dec( y3)  we  calculate  the  estimates  for  the  gain 
parameters,  J32 ,  and  J3l3.  Utilizing  these  estimations,  we  make  a  new  post  calculation 
based  on  (1.7),  yielding 

s3  dec(y3  ^>,3^2  /^i.3^i)-  (1-12) 

For  the  four-receiver  system  where  we  define  s4  =  dec( y4) ,  the  gain  estimates  are 
/?3  4 ,  /32A  and  /?,  4 .  These  estimations  are  now  used  in  the  post  calculation,  or 

s4  =  dec( y4  -J33A s3  -j32As2  -j3lA sx).  (1.13) 

With  the  formulas  ready  for  use,  the  Monte  Carlo  simulations  involve  first 
iterating  through  different  gain  parameter  values  with  fixed  SNR  levels.  The  estimates 
are  compared  to  the  actual  values  of  a  .  These  simulations  are  also  conducted  at  various 
SNRs  to  observe  the  effects  of  noise  on  the  estimators.  From  the  simulation  results,  we 
observe  that  as  any  of  the  a  gain  parameters  approach  1.0,  the  performance  of  the 
estimators  begins  to  degrade.  The  performance  is  degraded  as  additional  receiver  noise  is 
added  to  the  system.  In  our  simulations,  we  also  investigate  fixing  the  a  parameters 
while  the  SNR  is  varied  to  different  levels.  Improving  accuracies  of  the  estimates  are 
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observed  when  the  SNR  is  increased.  With  the  generation  of  the  post  calculation  signals, 
the  performance  of  the  RSIC  technique  is  illustrated  by  comparing  the  SER  between  the 

s  signals  and  theoretical  values  for  QPSK  modulation. 

It  is  made  clear  that  in  the  four-receiver  system  the  compounding  effect  of 
additional  interferences  has  a  negative  effect  on  estimator  performance.  Estimate 
accuracy  is  seen  to  decrease  with  the  slightest  addition  of  interference  amplitude  gain.  A 
method  to  mitigate  this  is  to  use  the  post  calculation  obtained  from  (1.13)  and  substitute 
those  values  into  a  new  calculation  of  the  gain  estimate.  New  (3  parameters  are  obtained 

and  can  again  be  used  to  calculate  a  new  s4 .  These  additional  rounds  of  estimate  re¬ 
calculation  are  seen  to  greatly  improve  the  accuracy  of  the  /?  parameters  as  well  as  SER 
performance  of  the  new  s4  signal. 

A  useful  metric  in  each  of  these  systems  is  the  signal-to-interference  ratio  (SIR) 
that  is  present  due  to  the  extraneous  transmissions  from  the  other  BSs.  It  is  useful  to 
compare  the  accuracy  of  the  estimates  as  a  function  of  SIR  to  see  how  overall 
interference  in  the  system  affects  the  performance  of  the  estimators.  Specifically,  in  the 
four-receiver  system,  it  is  clear  that  SIR  plays  a  large  role  since  there  are  three  signals 
interfering  with  the  SOI.  Even  small  values  of  amplitude  gain  from  one  of  these  signals 
make  a  considerable  effect  on  performance  when  compounded  with  the  gains  from  the 
others.  It  is  shown  in  the  simulation  results  that  as  the  SIR  for  a  specific  system 
increased,  the  estimator  performance  improves  and  in  turn,  the  SER  of  the  post 
calculations  is  diminished. 

RSIC  is  an  effective  method  for  extracting  a  desired  signal(s)  from  a  receiver(s)  in 
an  interference-filled  environment,  provided  of  course  that  the  estimates  calculated  for 
the  amplitude  gain  of  interference  signals  are  accurate.  Favorable  placement  of  each 
receiver  is  essential  to  minimize  the  effect  of  interference  gain  and  maximize  the  received 
strength  of  the  SOI,  so  that  the  estimators  in  the  system  perform  effectively.  This  thesis 
successfully  provides  a  method  of  providing  accurate  estimates  in  multiple  systems 
located  in  interference-filled  environments. 
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I.  INTRODUCTION 


A.  OVERVIEW 

The  increased  use  of  cellular  technology  has  provided  an  unparalleled 
advancement  of  communication  in  today’s  world.  Cell  phone  towers  and  base  stations 
(BS)  number  in  the  hundreds  in  just  a  small  geographic  area  making  signal  interference 
an  unavoidable  consequence.  While  several  methods  are  currently  used  to  mitigate  the 
amount  of  attenuation  and  multipath  interference  to  provide  clearer  links,  another  method 
exists  that  allows  the  removal  of  interferences  from  a  received  signal  so  that  only  the 
desired  information  remains.  This  method  is  accomplished  with  the  use  of  multiple 
receivers  positioned  so  that  an  initial  reference  signal  is  collected  that  is  passed  on  to  the 
next  receiver  to  effectively  conduct  interference  cancellation  (IC).  Consequently,  this 
method  does  not  allow  an  immediate  demodulation  of  the  clean  signal.  In  a  multi¬ 
platform  scenario,  with  more  receivers  involved  in  the  cancellation  procedure,  more 
latency  and  error  is  expected  in  the  final  signal  demodulation.  Even  so,  the  additional 
time  and  processing  power  necessary  to  execute  this  procedure  will  still  benefit  those 
conducting  long-term  surveillance  and  collection,  where  some  latency  is  acceptable  or 
real-time  signal  processing  may  not  be  required. 

Successive  interference  cancellation  (SIC)  is  considered  a  member  of  the  family 
of  post-IC  techniques  which  are  focused  on  the  cancellation  or  suppression  of 
interference  on  an  a-posteriori  basis  [1].  This  means  that  the  low  complexity  SIC  method 
must  first  collect  a  signal  in  order  to  then  ’’demodulate  and/or  decode  desired 
information,  and  use  this  information  along  with  channel  estimates  to  cancel  received 
interference  from  the  received  signal”  [2].  A  related  SIC  scheme  known  as  reference- 
based  successive  interference  cancellation  (RSIC)  makes  use  of  a  multi-platform  system 
and  what  is  called  a  reference  signal  to  eliminate  interference  in  received  signals. 

Previous  work  concerning  RSIC  [3]  explains  how  by  placing  a  receiver  in  a 
favorable  location,  a  signal  is  collected  and  demodulated  so  that  it  is  mostly  interference- 
free  (now  known  as  a  reference  signal)  and  then  transmitted  to  the  following  receiver  in 
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this  multi-platform  system.  This  next  receiver,  which  is  positioned  so  that  it  collects  a 
second  signal  interfered  with  a  scaled  version  of  the  first  signal  (scaled  due  to  amplitude 
gain  or  loss),  can  now  use  that  clean  reference  signal  and  a  calculated  estimate  of  the  gain 
to  subtract  the  interference  from  its  overall  received  signal.  Now  a  cleaned  up  version  of 
the  second  signal  is  available  and  is  also  called  a  reference,  but  for  the  second  signal. 
This  reference  signal,  in  addition  to  the  one  from  the  first  receiver  is  passed  on  to  a  third 
receiver.  The  cancellation  procedure  is  repeated  and  a  cleaned  up  reference  for  signal 
three  is  available  to  a  fourth  receiver. 

Thus,  the  procedure  applies  to  multiple  receivers.  In  practice,  the  subsequent 
receivers  are  the  ones  to  calculate  the  amplitude  gain  estimates  for  any  interfering  signal 
and  to  conduct  the  cancellation  to  allow  the  signal  of  interest  (SOI)  demodulation.  To 
note,  the  study  in  [3]  simply  utilizes  various  arbitrary  values  for  the  amplitude  gains  and 
estimates  solely  to  simulate  and  observe  the  performance  of  the  cancellation  technique. 
The  purpose  of  this  thesis  is  therefore  to  model  several  communication  systems  with 
multiple  interference  configurations  and  estimate  the  gains  of  the  received  interfering 
signals  and  evaluate  the  performance  of  the  RSIC  technique  in  terms  of  symbol  error  rate 
(SER)  of  the  SOI  using  quadrature  phase  shift  keying  (QPSK)  modulation.  In  addition, 
we  evaluate  the  accuracy  of  the  estimates. 

B.  SCOPE  AND  ORGANIZATION 

In  order  to  successfully  estimate  these  gains,  a  mathematical  signal  model  is 
needed  to  describe  the  received  signals  and  interferences  in  the  multiple  systems.  Next, 
the  least  squares  error  (LSE)  method  of  estimation  is  performed  to  calculate  the 
amplitude  gain  parameters  of  the  unwanted  signals.  These  calculations  are  compared  to 
the  actual  gain  values  from  the  interference  and  then  the  reference  signals  are  computed 
with  the  help  of  these  estimates  and  passed  on  to  the  next  receiver  in  order  to  continue  the 
RSIC  technique.  Utilizing  Monte  Carlo  simulations,  we  calculate  SER  against  various 
signal-to-noise  ratios  (SNR). 

In  Chapter  II,  the  concept  of  how  the  multi-platform  system  uses  the  RSIC 
technique  is  discussed  in  detail.  The  necessary  gain  parameter  estimators  are  calculated 
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and  the  accuracy  of  the  estimates  is  investigated  for  a  two-receiver  system  using  QPSK 
modulation.  The  results  of  Monte  Carlo  simulation  for  a  two-receiver  system  are 
examined  by  observing  how  the  estimates  change  while  first  holding  different  SNRs 
constant  followed  by  next  holding  the  gain  parameters  constant. 

In  Chapter  III,  a  third  receiver  is  added  to  the  system  and  the  mathematical 
formulas  needed  for  the  parameter  estimation  are  discussed.  The  performances  of  the 
estimators  for  a  three-receiver  QPSK  system  are  reviewed. 

In  Chapter  IV,  a  fourth  receiver  is  added  and  the  calculations  for  the  four-receiver 
system  are  discussed  and  the  system  is  simulated.  The  performances  of  the  estimators  for 
a  four-receiver  QPSK  system  are  then  discussed. 

In  Chapter  V,  we  present  our  conclusions  and  provide  recommendations  for 
follow  on  work. 
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II.  INTERFERENCE  CANCELLATION  AND  PARAMETER 
ESTIMATION  FOR  A  TWO-RECEIVER  SYSTEM 


A.  REFERENCE-BASED  SUCCESSIVE  INTERFERENCE  CANCELLATION 

RSIC,  or  SIC  for  short,  is  a  technique  that  removes  unwanted  signal  interferences 
from  a  set  of  receivers  in  order  to  obtain  one  or  multiple  SOIs  from  a  set  of  transmitters. 
These  transmitters  can  exist  in  one  of  several  configurations.  A  good  example  of  a 
configuration  is  shown  in  Figure  1,  where  the  coverage  areas  of  multiple  cellular  BSs 
overlap  with  those  from  subsequent  BSs,  resulting  in  a  layout  where  different  regions 
contain  multiple  signals.  Certain  sectors  of  the  coverage  areas  exist  that  contain  only  a 
specific  number  of  these  signals. 


Base  Station  Sectoring  Coverage  Area 


Figure  1.  Configuration  of  multiple  BSs  with  overlapping  antenna  coverage  which  has 
different  signal  combinations  received  in  each  sector.  From  [3]. 

Each  sector  is  modeled  so  that  it  contains  a  received  signal  y;  that  is  described 


by: 


y,  =5,+^, 


(2.1) 


(2.2) 


_y3  —  s3+a2  3s2  +  ai  35'1  +  w3 , 


(2.3) 


y 4  S4  ~\~  CC 3  4^3  +  ^2,4^2  ^1,4^1  ^4  ’ 
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(2.4) 


where  amn  is  the  amplitude  gain  (or  loss)  of  the  signal  from  BS  m  received  by  receiver  n 

and  is  assumed  real  and  positive,  and  where  w  is  the  additive  white  Gaussian  noise 
(AWGN)  in  the  receiver.  For  purposes  of  this  thesis,  the  energies  of  each  QPSK  signal  y 
for  i  =1,2, 3, 4  are  made  identical,  the  signals  are  assumed  to  have  perfect  phase 
synchronization  with  one  another  and  the  a  values  define  their  amplitude  ratio. 
Receivers  in  each  sector  are  assumed  to  only  have  the  knowledge  of  the  signal  y  they 
collect  and  the  reference  signals  sent  by  the  other  receivers  (i.e.,  receiver  two  (Rx2)  only 
knows  y2  and  y  ;  Rx3  only  knows  y3 ,  s2  and  y ;  Rx4  only  knows  y4 ,  s3 ,  s2  and  v, ).  If 

for  example  y  is  the  SOI  transmitted  by  base  station  four  (BS4),  then  Rx4  placed  in 
Sector  D  picks  up  the  SOI  as  well  as  interference  from  the  three  other  BSs  (i.e.,  y4  is  the 

received  signal  in  Sector  D).  Analytically,  each  one  of  these  received  signals  has  a  signal- 
to-interference  ratio  (SIR)  that  affects  the  performance  of  the  cancellation  technique. 
(This  concept  is  discussed  in  detail  in  Section  C  of  this  chapter.)  SIC  is  implemented  with 
Rxl  located  in  Sector  A  collecting  an  interference-free  or  clean  version  of  y  and 
transmitting  it  to  Rx2  located  in  Sector  B.  Next,  Rx2  conducts  the  calculations  needed  for 
gain  parameter  estimation  and  subtracts  or  cancels  the  scaled  version  of  y  from  the 

received  signal  y2 ,  resulting  in  an  estimate  for  y  (a  new  reference  signal).  Rx2  now 
transmits  this  new  reference  signal  in  addition  to  the  first  reference  signal  over  to  Rx3. 
This  iteration  of  cancellation  continues  until  only  the  clean  estimate  of  y  is  left. 

Clearly  for  this  method  to  effectively  work,  the  receivers  need  to  have  a  good 
estimate  of  the  values  of  ai  .  to  properly  conduct  the  IC.  This  gain  parameter  estimate  is 
called  /3  ,  which  is  eventually  used  for  the  estimated  signals  given  by 


=  Ai  —  ^ , 

(2.5) 

■y  =  y2~P\,2s  i  ~w2> 

(2.6) 

S3  y3  P2,3S2  Pl,3Sl  W3’ 

(2.7) 

•A  =  T4  -  A,4'V3  -  A,4'S2  -  A,4'Si  -  W4, 

(2.8) 
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where  w  is  the  noise  in  the  receiver.  For  high  SNR,  the  noise  is  considered  negligible  and 
thus  the  receiver  becomes  interference-limited.  For  low  SNR,  the  receiver  becomes 
noise-limited  and  therefore  the  noise  may  even  have  a  greater  effect  on  receiver 
performance  than  the  interferences  themselves.  In  our  simulation  results  (presented  in 
each  chapter),  varying  cases  of  noise  and  interference  powers  are  considered.  The 
previous  work  on  this  topic  [3]  explored  the  performance  of  various  types  of  IC,  the  best 
of  which  is  demodulated  RSIC  where  the  receiver  first  demodulates,  or  makes  the 
decision,  as  to  what  information  is  collected.  This  clean  signal  is  then  passed  on  to  the 
next  receiver  to  continue  the  SIC  process.  The  decision-making  process  is  discussed  in 
the  next  section. 

B.  PARAMETER  ESTIMATION  FOR  A  TWO-RECEIVER  SYSTEM 
1.  Least  Squares  Error  (LSE)  Method 

Starting  with  the  notion  of  a  simple  two-receiver  system  using  Equations  (2.1) 
and  (2.2),  we  utilize  the  LSE  approach  from  [4]  to  determine  an  estimate  for  the  a 
values.  In  this  thesis  we  assume  discrete-time  signal  models.  This  is  appropriate  as 
receiver  signals  are  almost  always  sampled  prior  to  signal  processing.  The  equation  to 
determine  the  estimation  of  the  amplitude  gain  first  seeks  to  minimize  the  magnitude  of 
squared  difference  between  the  received  signal  y2  and  the  combined  data  of  s2  +  ax  2.v, .  It 

is  given  by 

JV-1  2 

J{ct  12)  =  -  02M  +  aX2sx[n])\  (2.9) 

n= 0 

where  J (a)  is  clearly  the  square-error  magnitude,  v2 [/zj,  y[n]  and  ,s2[/zj  are  the  discrete¬ 
time  complex  signals  collected  by  Rx2  which  is  located  in  Sector  B.  If  we  let 
y2  =[y2[0]y2[l]...y2[7V-l]]r  and  s,.  =  [y[0]y  [l]...y[N-ljf  for  i=  1,2  we  get 

J(ah2 )  =  |y2  -  (s2 + «i,2si  )f  =  (y2  -  s2  -  «i,2si )"  (y  2  -  s2  -  «usi ).  (2.10) 

where  (-)H  is  the  Hermitian,  or  complex  conjugate  transpose  operation.  Expanding 
Equation  (2.10)  yields 
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(2.11) 


J(a- i,2)  =  y2"y2  -y2"s2  -«i,2y2fls.  -s2"y2  +s2fls2  +al2s2"s1 

~  ^l,2Sl  V2+«i,2S|"s2+«i22s"S| 


and  taking  the  partial  derivative  of  (2.1 1)  with  respect  to  al  2  yields 


dJ 

da{: 


y2  Sj+s2  s,  Sj  y2+s,  s2  +  2cTj  2Sj  sr 


Setting  the  derivative  to  zero  and  then  solving  for  the  estimate,  we  get 

2^1, 2si  Si  =y2  si  +si  y2-s,  s2  -s2  Sj 


(2.12) 


(2.13) 


and  thus 


a 


1,2 


y/^+s^y.-s/^-s/s, 

2sl//s1 


(2.14) 


Using  the  definition  of  a  real  number,  we  can  reduce  Equation  (2.14)  to 

2Re{y2//s1}-2Re{s1//s2} 


a. 


u  2SjffS| 


(2.15) 


or  simply 


<xxl 


Re{y2%}-Re{A2} 


(2.16) 


where  px  l  =s,"s2  and£Si 

The  first  hurdle  in  this  estimation  problem  is  apparent  in  (2.16).  In  order  to  solve 
for  al2,  we  need  Sj  and  s2  which  are  currently  unknown.  The  best  chance  of  obtaining 

the  first  clean  reference  signal  is  with  the  receiver  placed  in  Sector  A.  Here,  Rxl  only  has 
its  own  internal  noise  and  chances  of  a  successful  demodulation  are  reasonable  given  a 
high  SNR.  Once  received,  the  signal  y:  is  sent  through  a  maximum  likelihood  detector 
(MLD)  [5]  in  Rxl,  which  estimates,  or  decides ,  which  symbols  are  obtained.  This 
estimate  made  from  (2.5)  is  designated  as  s, ,  that  is 
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Sj  =  dec(  yj)  =  dec{sx  +  wx). 


(2.17) 


where  dec{ )  signifies  MLD  operation.  Next,  the  decision  as  to  what  symbols  are  received 
by  Rx2  is  made  based  on  y2  alone,  such  that  we  have  an  initial  (pre)  calculation  for  s2, 

or  simply  s2  where 

s2=dec(y2).  (2.18) 

This  estimate  is  then  used  in  (2.16)  to  determine  an  estimate  for  a]2,  which  is  given  as 

_  Rely^s.J-Rejs.^s,,} 

fl.2  =  2  l:H,  1  2  (2-19) 

Si  s. 

The  value  for  J3l  2  is  now  used  to  compute  a  final  (post)  calculation  for  the  QPSK 
symbols  of  s,  in  (2.6),  or 

s2  =  dec( y2  -/?I2s,).  (2.20) 

2.  Cramer-Rao  Lower  Bound  (CRLB) 

The  CRLB  is  calculated  in  order  obtain  a  lower  bound  on  the  variance  of  our 
estimator  in  (2.19).  This  is  to  determine  whether  or  not  the  estimator  is  the  mean  value 
unbiased  estimator,  which  would  indicate  that  (2.19)  has  the  lowest  variance  of  any  other 
unbiased  estimator  for  all  possible  values  of  al2.  Starting  with  a  general  approach  from 

[4]  of  a  real  signal  in  white  Gaussian  noise,  we  assume  that  a  deterministic  signal  with  its 
dependence  on  an  unknown  parameter  0  is  observed  in  AWGN  as 

y\n\  =  s[rv,  0\  +  w[n]  n  =  0,1,. ..,iV  —  1.  (2.21) 

The  probability  density  function  (PDF)  is 

p(y;  0)  =  — ■  jv/  exp  j-  X(yW  -  °^2  \ ■  (2-22) 

(2;rcr2)22  l  2cr  «=o  J 

Taking  the  expected  value  of  the  second  derivative  of  the  log-likelihood  function,  we  get 
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J  d2  In  p(y;  0))  =  1  ^fds[rv,G] 

{  SO2  )  a2  til  dO 


(2.23) 


so  that  the  CRLB  on  the  variance  of  the  estimator  9  of  6  is  expressed  as 


var($)  > - ,  0  ^  = - — - v  • 

J  d~  In  p(y;0)  ^  ds[rv,Ot 

l  d°2  J  til  d0  , 


(2.24) 


For  use  in  this  thesis,  if  the  complex  terms  from  (2.2)  are  used  and  we  allow 
,v[  it,  0]  =  ax  2s\n\ ,  then  the  PDF  for  complex  valued  signals  is 


P( y2;«i,2)  =  ,exp{-(y2 -s2-al2s,)wC;'(y2 -s2  -a^s,)}.  (2.25) 


If  we  assume  that  Cy  =  a2 1  andCy  1  -  —  /  ,  the  log-likelihood  yields 


f  \  \  i 

lnp(y2;«i2)  =  ln  +  In- exp  — 2(y2-s2-a12s1)H(y2-s2-«12s1) 

\n  <7  J  [  a~ 

(  l  3  i 

=  ln  —  in  — 2(.y2Hy2-y2Hs2-al2y2Hsl-s2Hy2+s2Hs2 

^7T  (7  y  (7 


+  al2 s2  Sj  -«12Sj  y2  +  «lj2s1  s2  +  a,i2Sj  s,). 


(2.26) 


We  assume  that  s,  and  s2  are  known  at  this  point.  In  practice,  these  are  merely  estimates 
themselves  and  contribute  to  the  overall  variance  of  the  estimate.  Taking  the  first 
derivative  with  respect  to  ax  2  gives 


alnfY;0V2)  = - 4(-y  A  -s,"y2  +s,"s2  +2«us1"Sl), 

oax  2  cr  v  ' 


(2.27) 


and  a  second  differentiation  results  in 


<32  In  p(y2',al2) 


-Y(s,"s,). 


(2.28) 
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The  expected  value  of  this  provides 


^ d2  in  p(y2;al2)^ 


da. 


1,2 


2-  s, 


(2.29) 


which  yields  a  CRLB  of 

var(al2)> 


1 


f  a  2 


-E 


d  In p(y2,a12) 

dal  2 


2-  Sj 


(2.30) 


where  cr2  is  the  variance  caused  by  the  random  AWGN  and  al2=(3{2.  Therefore, 
Equation  (2.30)  shows  the  variance  of  /?,  2  is  in  terms  of  the  energy  of  the  first  signal  and 
is  effectively  the  theoretical  value  as  if  there  is  no  s2  present. 


C.  RESULTS  FOR  A  TWO  RECEIVER  SYSTEM 
1.  Fix  S2NR,  Vary  a12 

As  y2  is  made  up  of  signal  s2  and  the  interference  of  the  scaled  form  ofs, ,  the 
SIR  is  therefore  defined  as 

E  II  II2 

SIR  =  — —  =  ™  .  (2.31) 

F  II  II2 
Pi  2^1 

With  the  energy  of  each  signal  normalized  to  have  the  same  value,  the  equation  that 
defines  the  relationship  between  s2  and  the  interference  from  Sj  simplifies  to 

S2IR  = -4-.  (2.32) 

«1,2 

Furthermore,  the  SNR  of  each  signal  is  varied  by  simply  changing  the  energy  of 
the  noise  power  Pn  in  the  received  signals.  In  our  Monte  Carlo  simulation,  we  generate 
two  random  QPSK  symbol  streams  combined  with  accompanying  noise  and  interference 
to  create  (2.1)  and  (2.2).  The  values  of  at ,  that  are  examined  range  from  0.01  to  1.0  at 
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steps  of  0.01  while  setting  the  SNR  for  s2,  or  S2NR,  at  multiple  levels  to  observe  the 
effect  on  the  accuracy  of  the  gain  estimate  px  , .  With  their  inverse  relationship,  as  2  is 
increased,  the  value  of  S2IR  decreases.  Initially,  Pn  is  set  very  low  in  order  to  first  model 
a  system  where  there  is  a  very  high  S2NR,  which  is  then  followed  by  modifying  Pn  to 
have  values  closer  to  unity  (i.e.,  equal  to  normalized  signal  energy)  to  model  low  S2NR. 
From  Figure  2,  we  see  that  the  dashed  blue  line  represents  the  actual  at  2  amplitude  gain 

parameter  used  in  the  simulation.  At  high  S2NR  levels,  the  calculated  estimate  values  lie 
closely  along  the  theoretical  line  until  al2  approaches  1.0.  If  the  noise  is  increased  and 

the  S2NRS  are  changed  to  low  levels,  the  performance  of  the  estimator  begins  to  decline. 
Even  at  minimal  values  of  al2,  the  estimates  deviate  from  the  line  representing  no 

interference  and  the  deviation  becomes  significant  when  Pn  =0.8.  Interestingly,  as  these 
estimates  get  worse,  they  begin  to  converge  towards  the  value  of  0.5.  This  is  most  likely 
due  to  the  fact  that  the  estimator  from  (2.19)  is  more  accurate  for  systems  where  E  is 

much  greater  than  E  ,  and  therefore  does  not  yield  good  estimates  for  high  «,  2  levels. 

Another  way  to  see  how  effective  the  estimates  are  for  /?,  2  is  to  plot  them  against 
the  actual  al2  parameter  values  as  displayed  in  Figure  3.  Again,  the  dashed  blue  line 
represents  actual  amplitude  gain  parameter  values  and  the  subsequent  curves  show  the 
estimates  at  various  S2NRS.  At  low  Pn  levels,  the  estimates  are  fairly  accurate  until  aI2 

reaches  a  value  of  about  0.7.  At  higher  cr12  values,  the  simulation  results  show  estimates 

that  again  diverge  from  the  theoretical  line  and  approach  a  value  of  0.5.  Thus,  we  can 
appreciate  the  effect  of  increasing  noise  and  interference  on  the  accuracy  of  the  amplitude 
gain  estimations. 
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(a)  (b) 

Figure  2.  Curves  of  f}x  2  versus  (vs.)  S2IR  at  various  noise  levels:  (a)  Performance 
curves  with  high  S2NR;  (b)  Performance  curves  with  low  S2NR. 


Figure  3.  Curves  of  J3l2  vs.  al2  at  various  noise  levels:  (a)  Performance  curves  with 
high  S2NR;  (b)  Performance  curves  with  low  S2NR. 
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2.  Fix «12,  Vary  S2NR 

A  second  approach  to  observing  the  performance  of  the  estimators  is  to  hold  the 
interference  gain  values  constant,  and  while  the  SNR  of  s2  is  increased,  examine  the 
performance  of  the  SIC  technique.  The  simulation  is  conducted  at  different  al2  values 

while  decreasing  the  amount  of  noise  in  the  receivers  and  thus  we  are  able  to  calculate  the 
SER  for  the  increasing  S2NR.  It  is  important  to  remember  that  on  the  figures,  the  pre 
calculations  are  the  estimates  made  from  making  a  decision  based  solely  on  y2 ,  or  (2.18), 
whereas  the  post  calculations  from  (2.20)  provide  better  signal  estimates.  From  Figure  4, 
it  is  seen  that  at  low  al2,  both  the  SER  for  the  pre  and  post  calculations  of  the  s2 

symbols  closely  follow  the  theoretical  values  obtained  for  QPSK.  Also,  due  to  the 
relationship  from  (2.32),  a  better  performance  is  therefore  expected  for  higher  S2IR 
levels.  When  the  amount  of  interference  is  increased  (i.e.,  increased  amplitude  gain  on 
the  energy  of  s, ),  the  performance  decreases  (i.e.,  SER  starts  to  suffer).  Clearly  the 
incorporation  of  J3{  2  in  the  IC  makes  an  impact  to  the  quality  of  the  cleaned  up  signal,  as 

pre  calculations  are  severely  affected  by  increased  interference  while  post  calculations 
show  only  a  slight  to  moderate  deviation  from  the  theoretical  line. 
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Figure  4.  Performance  curves  of  SER  vs.  S2NR  for  pre/post  decisions  of  s2  at  different 
interference  combinations:  (a)  al2  =  0.2,  S2IR  =  14  dB;  (b)  al2  =  0.4,  S2IR  =  7.96 
dB;  (c)  al  2  =  0.6,  S2IR  =  4.4  dB;  (d)  al2  =  0.8,  S2IR  =  1.94  dB. 


Using  the  same  approach  of  setting  different  al2  values  and  varying  SNR,  we  see  the 
accuracy  of  the  J3l2  estimates  shown  in  Figure  5.  By  plotting  a  line  of  the  actual  gain 
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value  at  2  and  then  showing  how  the  estimates  approach  that  line,  a  clear  understanding 

is  made  of  how  greater  S2NR  plays  a  role  in  parameter  estimation.  Estimates  obviously 
improve  with  reduction  of  noise  in  the  system.  However,  with  higher  interference  levels 
such  as  when  ax  2  =  0.8,  increased  S2NR  does  not  provide  the  same  level  of  accuracy 

previously  seen  for  the  gain  estimate  at  lower  al2  values. 


(c>  (d) 

Figure  5.  Plots  of  the  actual  gain  parameter  ax  2  and  the  estimator  fil2  vs.  S2NR  values 
for  multiple  simulations  of  y2 :  (a)  a12  =  0.2,  S2IR  =  14  dB;  (b)  a12  =  0.4,  S2IR  = 
7.9  dB;  (c)  a12  =  0.6,  S2IR  =  4.4  dB;  (d)  al  2  =  0.8,  S2IR  =  1.9  dB. 


If  we  refer  back  to  the  CRLB,  it  is  possible  to  plot  the  variance  of  the  fil2 

estimates  calculated  from  (2.30)  against  the  theoretical  variance  as  if  the  system  is 
without  interference.  From  Figure  6,  it  is  seen  that  while  ax  2  remains  low,  it  is  possible  to 
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Variance  Variance 


attain  a  very  low  variance,  which  corresponds  to  a  more  accurate  estimate  2 .  However, 
once  aX2  starts  to  approach  1.0,  even  very  high  levels  of  SiNR  result  in  poor  variance  of 
the  estimate. 


Figure  6.  Plots  of  the  variance  of  the  estimator  /3X  2  vs.  SiNR  values:  (a)  ax  2  =  0.2,  S2IR 
=  14  dB;  (b)  al  2  =  0.4,  S2IR  =  7.9  dB;  (c)  a]  2  =  0.6,  S2IR  =  4.4  dB;  (d)  al  2  =  0.8, 

S2IR=  1.9  dB. 


With  the  basic  understanding  of  how  RSIC  is  applied  and  performs  for  two 
receivers,  we  next  observe  the  IC  technique  for  the  simulation  of  a  three-receiver  system. 
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III.  ADDITION  OF  A  THIRD  RECEIVER 


Continuing  with  the  idea  of  successive  receivers  within  a  multi-platform  system, 
we  assume  that  a  third  signal  is  obtained  by  Rx3  located  in  Sector  C  and  contains 
information  from  BS1,  BS2  and  BS3,  as  seen  in  Figure  1.  Clearly  the  procedure  involved 
for  this  estimation  is  more  involved  as  the  received  signal  y3  contains  two  interfering 
signals. 


A. 


PARAMETER  ESTIMATION  FOR  A  THREE-RECEIVER  SYSTEM 


Following  the  same  method  from  Chapter  II,  we  see  that  a  third  receiver  would 
make  use  of  (2.3)  and  the  LSE  approach  gives 


J(a)  =  y,  y3[n\-(s3[n\  + a23s2[i'i\  + ax3sx[n\)\  . 


n= 0 


(3.1) 


Taking  the  vector  form  of  (3.1)  and  expanding,  we  get 


J(a)=  y3 -(s3 +a2j3s2 -t-ar^Sj) 


=(y3 


■  3s2  —  otj  3S1 )  (y3  —  s3  —  oc~,  3s2  —  ctj  3Sj ) 


(3.2) 


and 


J (cc')  y3  y3  y3  s3  0^2, 3y 3  s2  ^c3y 3  Sj  s3  y3~^s3  s3  +  cr23s3  ^2 


H 


2  „  H„ 


+  CC\  3S3  Sj  — ^2  3^2  y3"^^2  3^2  S3  +  (X^  3^2  ^2  ^"1  3^2  3^2  ® 
—  Ct\  3S1  y^+^^Sj  S3  +  GC\  36^2  3S1  S2  H"  CL\  3Sj  Sj. 


(3.3) 


The  important  difference  from  the  two-receiver  system  is  the  need  to  now  estimate  for 
two  gain  parameters.  Taking  the  partial  derivatives  with  respect  to  the  individual  gains, 
we  get 


dJ 

da23 


y3  s2 


s2  s2  y3 ■ 


■  2«23s2 


dJ 

dax  3 


-y,"s 


s  ws 
1  1  a3  ai 


^2,3^2 


-N  y3  + 


Sf  s3 


■  «2,3S|  S2  ■ 


■  2 ax  3SjHsr 


(3.4) 


19 


Setting  to  zero  and  solving  for  each  one  of  the  estimates  we  get 


a2.3 


y3  s2+s2  y3  s2  s2  s3  cqj3(s2  Sj +Sj  s2) 
2s2hs2 


a13  — 


y3"s,  +  s"y3  -Sj^Sj  -s"s3  -a23(s2HSj  +8,%) 


(3.5) 


2Si//s1 


We  can  utilize  the  definition  of  a  real  number  to  simplify  the  equations  and  get 


2  Re  {y  3Hs2 }  -  2  Re{s,ffs2 }  - 


^2,3 


2Re{y,  Si}-2Re{s3  Si}-a23(2Re{s2  s, }) 


2sihs1 


(2Re{s2\}) 


2s2  s2 


2Re{y3Hs2}-2Re{p32}- 


2Re{p21}Re{y3gs,}-2Re{>c»21}Re{p31}-a23(2Re{p21}2) 

E. 


(3.6) 


2  E. 


£s,Re{y3HS2}-£s,Re{p3,2}-Re{y92,i}Re{y3M  +  Re{P2,i}Re{P3.i}  +  «2.3Re{A 


«2.3  (^s,  ^ )  -  «2,3 Re ( /^2,1 } 2  =  ( )  ( Re  { y 3" s2 }  -  Re ( P3,2 } )  -  ( Re  ( /^2,1 } )  ( Re  { y 3" Si }  -  Re  { /?3 1 } ) . 


The  same  approach  is  taken  for  the  estimator  of  a{ ,  and  the  following  sequence  of 
equations  is  obtained: 
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2Re{y3\ }  -  2  Re{s3\ }  - 

2  Re{y3Hs, }  -  2  Re{s3Hs2 }  -  alj3(2  Re{s2Hs, } ) 

(2Re{s/Sl}) 

2s,Hs2 

2SjWSj 

2Re{y3ws1}-2Re{/731}- 

2  Re  { p2A }  Re  { y  " s 2 }  -  2  Re  { p2 s }  Re  { p3  2 }  -  at  3 

,(2Re{/72J}2)" 

E , 

2E \ 

£S2Re{y3M-£S2Re{p3ji 

}  -  Re  { p2 a }  Re  {y  3% }  +  Re  { p2  l }  Re  { p3  2 }  +  al3 

Re{/C2,,}2 

«l,3(£S|£S2)-«URe{/72.l}2=(£S2)(Re{y3Hsi}-Re{A.l})-(Re{/72.l})(Re{y3/,S2}-Re{A.2})- 

Now  solving  for  the  a  values  in  (3.6)  and  (3.7),  we  arrive  at  the  final  estimators  for  the 
amplitude  gains  of  a  three-receiver  system  which  are  given  by 

_  ESi  (Re{y3Hs2}-Re{/?32})-(Re{/?21})(Re{y3Hs1}-Re{/731}) 

^  £Si£S2-Re{A/  ’ 

(3.8) 

(^^-Relp,,}2)  ' 

There  is  however  the  same  problem  that  is  experienced  in  Equation  (2.16)  where 
the  calculations  in  (3.8)  require  the  signal  s3.  Therefore  the  system  must  first  make  a 

decision  as  to  what  symbols  are  transmitted  based  only  on  the  received  signal.  In  this 
scenario,  Rx3  makes  the  pre  decision  on  the  received  signal,  which  is  given  by 

s3  -dec(y3).  (3.9) 

Substituting  (3.9)  into  (3.8)  yields 

_  £s,  (Re{y3Hs2}-Re{A.2})-(Re{A.i})(Re{y3Hs1}-Re{A.i}) 

Pl3~  (EsES2-Re{p2J) 


(£sA-Re{/32,!}2) 
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(3.10) 


where  p32  =S3HS2,  p3X  =  S3HSX  and  p2l  =S2HSl.  Utilizing  the  calculated  gain  estimates, 
we  make  a  new  post  calculation  based  on  (2.7)  yielding 

s3  =  dec(y3 - /323s2 (3.11) 


1.  CRLB  Calculation 

If  we  tailor  the  complex  PDF  from  (2.25)  for  a  three-receiver  system,  we  get 
p( y3;a)  =  |exp{-(y3  -s3  -a2>3s2  -ax3sx)HC~l{y3  -s3-a23s2-ax3sxj\.  (3.12) 


The  second  derivative  of  the  log-likelihood  for  each  gain  parameter  yields 

<32  In  p(y3\a23)  2 


da 


2,3 


d~\np(y3;al3)  _  2  f  H  ^ 

"""  dal  ~  S'*' 


from  which  the  CRLB  for  each  parameter  is 


(3.13) 


var(  a2 , )  > 


1 


-E 


f  d2\np(ypa23) 

V  ^^2,3  J 


2-  s. 


var(«]  3)  > 


1 


f  ^2 


-E 


d~  In  p(y2',al  3) 


da[ 


3  J 


2-  s, 


(3.14) 


where  a2  3  =  fi23  and  aX3  =  fil  3 . 

As  can  be  seen  from  (3.14),  the  theoretical  variance  of  an  estimate  solely  depends 
on  the  noise  of  the  system  and  energy  of  the  signal  corresponding  to  the  amplitude  gain. 
The  actual  variance  calculated  in  the  simulation  will  depend  on  the  overall  interference 
experienced  at  the  receiver. 


22 


B. 


RESULTS  FOR  A  THREE-RECEIVER  SYSTEM 


In  order  to  best  illustrate  the  effect  of  different  amplitude  gains  of  the  interfering 
signals,  the  simulations  are  performed  with  the  assumption  that  Rx3  has  access  to  the 
good,  clean  reference  signals  s,  and  s2 .  Therefore,  in  the  simulations  y,  and  y2  from 
Equations  (2.1)  and  (2.2)  are  assumed  to  have  high  SNRs  and  SIRs  for  their  respective 
signals.  For  Equation  (2.2),  this  means  assuming  that  s,  includes  only  a  small 

interference  gain  so  that  when  a  decision  is  made  on  y2,  the  reference  signal  passed  to 

Rx3  is  mostly  error-free.  To  illustrate  the  result  of  the  SIC  technique  on  a  three -receiver 
system,  first  the  SNR  levels  are  fixed  and  the  cc  values  are  varied.  This  is  followed  by 
the  reverse,  through  fixing  the  a  values  and  varying  the  SNR  of  the  signals  and 
observing  the  performance  of  the  estimators.  With  estimation  needed  for  two  gain 
parameters,  the  simulation  is  set  up  so  that  while  the  SNR  of  the  system  is  fixed,  we 
iterate  through  several  different  interference  combinations  by  holding  one  gain  at  a 
certain  level  and  running  calculations  through  a  range  of  the  other  gain  values. 


1.  Fix  S3NR,  Vary  a2  3  and  a13 

With  two  signals  now  interfering  with  s3 ,  the  performance  of  the  estimators 

depends  on  the  signal-to-total-interference  ratio  of  the  system.  As  is  seen  with  (2.31),  this 
new  relation  is  defined  as 

=  i - iH - ¥•  <3'15) 

^s2  '  ^Sj  \\a23s2\\  +  3SJ 

where  s,  and  s,  are  considered  uncorrelated  based  on  being  random  QPSK  signals 
generated  by  two  different  transmitters. 

Since  the  energies  of  s, ,  s2  and  s3  are  the  same,  the  resulting  equation  yields 


SJR  = 


a23  +  ^1,3 


(3.16) 
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Unless  otherwise  denoted  as  an  individual  signal  interference  (i.e.,  S3FR  and  S3I1R), 
Equation  (3.16)  is  used  to  describe  the  total  signal-to-interference  ratio.  Consequently,  as 
S3IR  describes  the  overall  interference  in  the  system,  in  many  cases  there  is  no  need  to 
hold  al3  constant  and  then  vary  a23  since  the  resulting  S3IR  is  the  same  and  would 

generate  similar  performance  with  regard  to  SER  and  estimation.  For  example,  if  cr,  3  = 
0.5  and  a13  =  0.2,  the  overall  S3IR  is  l/(0.52  +0.22)  =5.4  dB  .  The  performances  of  the 
estimators  in  the  system  are  the  same  if  a23  =  0.2  and  ai3  =  0.5.  Further  examination  of 
this  relationship  is  seen  in  the  Appendix. 

The  first  way  to  gauge  the  effect  of  the  estimation  accuracy  is  with  various  noise 
levels  in  the  simulations  held  constant  and  having  the  gain  parameters  (and  therefore 
S3IR)  change  values.  Viewing  the  results  in  this  fashion,  we  observe  how  close  the 
calculations  of  (3.10)  are  compared  to  the  theoretical  values.  From  Figure  7,  it  is  seen  that 
as  a23  is  varied  at  high  S2NR  levels,  the  estimates  remain  accurate  with  low  al2  values 

such  as  0.01  and  0.2,  but  heavily  degrade  at  0.5  and  0.8.  From  Figure  8,  we  see  that  a 
three-receiver  system  suffers  greatly  in  environments  with  low  SiNR,  as  the  estimates 
barely  approach  the  theoretical  line.  If  the  estimates  are  plotted  against  the  actual 
interference  values,  a  similar  trend  is  observed  in  which  the  estimates  fall  back  to  the 
value  of  0.5  under  high  amounts  of  interference  as  seen  in  Figures  9  and  10.  From  these 
results,  we  see  that  both  noise  and  interference  maintain  a  profound  effect  on  the  outcome 
on  the  estimators  from  Equation  (3.10). 
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Figure  7.  Plots  of  gain  estimates  /?2  3  vs.  S3I2R  values  with  fixed  high  S3NR  levels: 
(a)  al3  =0.01;(b)  al3  =0.2;  (c)  al3  =0.5;(d)  a13  =0.8. 
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Figure  8.  Plots  of  gain  estimates  fi23  vs.  S3I2R  values  with  fixed  low  S3NR  levels: 
(a)  al3  =  0.01;  (b)  al3  =  0.2;  (c)  al3  =  0.5;  (d)  al3  =  0.8. 
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Figure  9.  Plots  of  gain  estimates  /?2  3  vs.  actual  a2  3  parameter  values  with  fixed  high 
S3NR  levels:  (a)  al3  =  0.01;  (b)  a13  =  0.2;  (c)  al3  =  0.5;  (d)  a13  =  0.8. 


(a) 


(t>) 


(C) 

Figure  10.  Plots  of  gain  estimates  J3-,3  vs.  actual  a23 
S3NR  levels:  (a)  «13  =  0.01;  (b)  al3  =  0.2; 


(«*i 

parameter  values  with  fixed  low 
(c)  aI3  =0.5;(d)  al3  =0.8. 
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2.  Fix  a, ,  and  a, 3,  Vary  S3NR 

To  observe  the  performance  differences  of  the  various  scenarios  in  a  different 
manner,  we  set  one  of  the  amplitude  gains  constant  and  the  second  is  varied  from  almost 
zero  to  moderate/heavy  interference.  The  simulation  is  performed  with  parameter  a2  3 

taking  the  values  0.01,  0.2,  0.5  and  0.8  while  the  parameter  a13  changes  to  0.01,  0.1,  0.2, 
0.4,  0.6  and  0.8  to  illustrate  their  respective  SERs  as  shown  in  Figures  11-14. 

Utilizing  the  same  pre  and  post  decision  concepts  from  Chapter  II,  we  see  the 
differing  performances  resulting  from  the  multiple  combinations  of  interference  gain 
values.  For  brevity,  only  the  post  SER  curves  are  discussed  since  these  are  the  results 
expected  after  the  full  application  of  the  SIC  technique.  Starting  at  a  low  a2  3 ,  we  see  that 

the  post  SERs  realized  with  high  S3NR  at  each  of  the  ax 3  values  are  quite  small,  around 
lO^-lO'6,  which  means  that  the  SIC  technique  works  well  and  provides  a  clean  s,  that  is 
close  to  the  original  signal.  This  occurs  until  the  interference  amplitude  from  the  first 
signal  is  set  to  0.8.  At  this  point,  the  SER  begins  to  degrade  as  seen  in  Figure  11.  If  the 
a23  is  increased  to  0.2,  then  the  performance  starts  to  deteriorate  earlier  at  al3  =  0.6  as 

shown  in  Figure  12.  From  Figure  13,  it  is  seen  that  SER  now  has  a  noticeable  decline  at 
al3  =  0.4,  and  from  Figure  14,  we  see  there  is  an  immediate  decrease  in  SER 

performance  at  every  S3NR  and  a13  level  due  to  the  high  value  of  a23  =  0.8.  At 
interference  levels  this  high,  it  is  very  difficult  to  successfully  retrieve  the  SOI. 
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Figure  1 1 .  Performance  curves  of  SER  vs.  S3NR  for  pre/post  decisions  of  s3  at  different 
interference  combinations  when  a23  =  0.01:  (a)  a13  =  0.01,  S3IR  =  37  dB;  (b)  al3 
=  0.1,  S3IR  =  20  dB;  (c)  «u  =  0.2,  S3IR  =  14  dB;  (d)  au  =  0.4,  S3IR  =  7.96  dB; 
(e)  a13  =  0.6,  S3IR  =  4.4  dB;  (f)  a13  =  0.8,  S3IR  =  1.94  dB. 
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(e)  (f) 


Figure  12.  Performance  curves  of  SER  vs.  S3NR  for  pre/post  decisions  of  s3  at  different 
interference  combinations  when  a2  3  =  0.2:  (a)  al3  =  0.01,  S3IR  =  14  dB;  (b)  a13 
=  0.1,  S3IR  =  13  dB;  (c)  al3  =  0.2,  S3IR  =  11  dB;  (d)  au  =  0.4,  S3IR  =  6.99  dB; 
(e)  a, 3  =  0.6,  S3IR  =  3.98  dB;  (f)  au  =  0.8,  S3IR  =  1.67  dB. 
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Figure  13.  Performance  curves  of  SER  vs.  S3NR  for  pre/post  decisions  of  s3  at  different 
interference  combinations  when  a23  =  0.5:  (a)  al3  =  0.01,  S3IR  =  6.02  dB;  (b)  al3 
=  0.1,  S3IR  =  5.85  dB;  (c)  a13  =  0.2,  S3IR  =  5.38  dB;  (d)  a13  =  0.4,  S3IR  =  3.87 
dB;  (e)  «l3  =  0.6,  S3IR  =  2.15  dB;  (f)  au3  =  0.8,  S3IR  =  0.51  dB. 
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Figure  14.  Performance  curves  of  SER  vs.  S3NR  for  pre/post  decisions  of  s3  at  different 
interference  combinations  when  a23-  0.8:  (a)  al3  =  0.01,  S3IR  =  1.94  dB;  (b) 
al3  =  0.1,  S3IR  =  1.87  dB;  (c)  au  =  0.2,  S3IR  =  1.67  dB;  (d)  «u  =  0.4,  S3IR  = 
0.97  dB;  (e)  «l  3  =  0.6,  S3IR  =  0  dB;  (f)  au  =  0.8,  S3IR  =  -1.07  dB. 
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Next,  the  accuracies  of  the  estimates  are  observed  alongside  the  actual  gain 
parameters  in  Figures  15-22.  Several  plots  at  varying  levels  of  interference  are  required 
to  show  how  changing  one  cc  parameter  affects  the  estimate  fi  of  the  other  parameter. 
In  every  figure,  each  subplot  indicates  a  different  run  of  the  simulation  with  gain 
parameters  fixed  to  certain  values.  For  example,  from  Figures  15  and  16,  the  initial  «,  3 

parameter  is  set  very  low  to  0.01  while  al3  moves  across  a  range  of  a  values.  The 
estimates  /?,  3  are  very  accurate  throughout  these  iterations,  however,  at  lower  S3I1R  the 
estimates  /?,  3  require  higher  S3NR  in  order  to  properly  calculate  the  gain  parameter.  The 
next  round  of  simulations  are  performed  with  a23  =  0.2  where  an  immediate  change  is 
observed  where  the  first  calculations  of  [5 , 3  lie  below  the  actual  value  line  as  seen  in 
Figure  17.  These  are  improved  with  a  higher  S3NR,  but  again  as  the  a13  interference 
approaches  1.0,  the  estimates  /?23  and  /3l3  begin  to  suffer  as  shown  in  Figure  18.  As 

both  the  interference  values  continue  to  intensify,  the  overall  accuracy  for  both  estimates 
decrease,  up  to  the  point  where  even  increasing  the  S3NR  actually  begins  to  have  a 
negative  effect  on  the  actual  estimation  values  as  shown  in  Figures  19-22. 

The  CRLB  for  a  three-receiver  shows  similar  results  that  are  obtained  in  the  two- 
receiver  system.  As  seen  in  Figures  23  and  24,  each  row  indicates  a  specific  configuration 
of  gain  parameters  for  the  simulation.  For  example,  the  top  row  shown  in  Figure  24  has 
a2  3  =  0.2  and  aX2  =0.1  whereas  the  bottom  row  is  setup  with  a2  3  =  0.2  and  a1  2  =  0.4. 

The  calculated  variances  from  each  of  the  simulations  follow  the  curve  of  the  theoretical 
variance  from  (3.14),  but  deviate  a  certain  degree  based  on  the  addition  of  interference 
into  the  system.  This  is  an  expected  result  as  additional  gain  from  interfering  signals 
effectively  degrades  the  performance  of  the  estimators,  causing  the  variance  to  increase 
in  each  estimation  calculation.  As  the  low  interference  simulations  from  Chapter  II  and 
III  confirm  the  theoretical  variance  values  of  (2.30)  as  well  as  of  (3.14),  the  CRLB  for  a 
four-receiver  system  will  not  be  calculated. 
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andpQ,  a,.,  and  aQ,andp 


a2  3  =  0.01 ,  S3I2R  =  40  dB  a,,  3  =  0.01 ,  S^R  =  40  dB 


(a) 


a2,3  =  0-01 ,  S3'2R  =  40  dB 


ai,3  =  0,1  ’  S3'iR=  20  dB 


(b) 


a23  =  0.01,S3l2R  =  40  dB 


ai,3  =  0-2,  S3'lR  =  14  dB 


(C) 


Figure  15.  Plots  of  the  actual  gain  parameters  a2  3 ,  cr13  and  the  estimators  3 ,  3  vs. 

different  S3NR  values  for  multiple  simulations  of  y3  while  a2  3  =  0.01:  (a)  S3IR  = 
37  dB;  (b)  S3IR  =  20  dB;  (c)  S3IR  =  14  dB. 
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and  a  andp  a,,andp 


0.2 


0  5  10  15 


S3NR  (dB) 


ai,3  =  0-4’  s3'iR  =  7-96  dB 


0  5  10  15 


S3NR  (dB) 


a  =  0.6,  S,I.R=4.44  dB 

1,3  3  1 


a2,3  =  001-S3l2R  =  40dB 


a1>3  =  °-8,  S3'iR  =  1  '94  dB 


Figure  16.  Plots  of  the  actual  gain  parameters  a2  3 ,  al3  and  the  estimators  f52  3 ,  /i,  3  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a2  3  =  0.01:  (a)  S3IR  = 
7.96  dB;  (b)  S3IR  =  4.44  dB;  (c)  S3IR  =  1.94  dB. 
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and  p0 ,  a^andp,,,  aQ.,andp 


a2,3  =  0-2>S3l2R=14dB 


^  3  =  0.01,83^  =  40  dB 


tt2,3  =  02’  S3'2R=14dB  a13  =  0.1,S3l1R  =  20dB 


a2  3  =  0.2,  S3I2R  =  14  dB  &1  3  =  0.2,  83!^  =  14  dB 


Figure  17.  Plots  of  the  actual  gain  parameters  a2  3 ,  at  3  and  the  estimators  f32  3,  ,  vs. 

different  S3NR  values  for  multiple  simulations  of  y3  while  a2  3  =  0.2:  (a)  S3IR  = 
14  dB;  (b)  S3IR  =  13  dB;  (c)  S3IR  =  11  dB. 
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3nd  P2.3  “2,3  aPd  ^2,3  “o  and  P 


“2,3  =  °'2’S3l2R=14dB 


“2,3  =  0-2’S3'2R=14dB 


“2,3  =  °'2’S3l2R=14dB 


“i,3  =  0-4’  ^R37-96  dB 


(a) 

a13  =  0.6,  83!^  =  4.44  dB 


(b) 

a,  3  =  0.8,  S^R-154  dB 


(c) 


Figure  18.  Plots  of  the  actual  gain  parameters  a2  3 ,  a13  and  the  estimators  /?23,  /?,  3  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a23  =  0.2:  (a)  S3IR  = 
6.99  dB;  (b)  S3IR  =  3.98  dB;  (c)  S3IR  =  1.67  dB. 
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andp0,  a,.,andp  a^andp 


a2  3  =  0.5,  S3I2R  =  6.02  dB 


0  5  10  15 


(a) 


S3NR  (dB) 


a2  3  =  0.5,  S3I2R  =  6.02  dB  ^  3  =  0.1 ,  S^R  =  20  dB 


a23  =  0.5,  S3I2R  =  6.02  dB 


ai,3  =  °'2,  S3>1R  =  14  dB 


Figure  19.  Plots  of  the  actual  gain  parameters  a2  3 ,  al3  and  the  estimators  ft2  3 ,  /?,  3  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a23  =  0.5:  (a)  S3IR  = 
6.02  dB;  (b)  S3IR  =  5.85  dB;  (c)  S3IR  =  5.38  dB. 
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andp0,  a^andp^  a,,andp 


a2  3  =  0.5,  S3I2R  =  6.02  dB 


ai,3  =  0-4’  s^R"7-96  dB 


a  =0.5,  S3I2R=6.02  dB 


a  =  0.6,  SJ  R  =  4.44dB 

1,3  3  1 


a23  =  0.5,  S3I2R  =  6.02  dB 


ai,3  =  0-8>  S3'iR  =  1  94  dB 


Figure  20.  Plots  of  the  actual  gain  parameters  a2  3 ,  a13  and  the  estimators  /?2  3 ,  fJl ,  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a23  =  0.5:  (a)  S3IR  = 
3.87  dB;  (b)  S3IR  =  2.15  dB;  (c)  S3IR  =  0.51  dB. 
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andp,,  a^andp^  a,,andp 


a23  =  0.8,S3l2R=1.94dB 


a13  =  0.01,S3l1R  =  40dB 


a23  =  0.8,S3l2R=1.94dB 


ai,3  =  O-1’  S3'iR  =  20  dB 


a2>3  =  °-8’S3l2R=1-94dB 


ai,3  =  °'2,  S3'lR=  14  dB 


Figure  21.  Plots  of  the  actual  gain  parameters  a2  3 ,  al3  and  the  estimators  /3-,3,  /?, ,  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a2  3  =  0.8:  (a)  S3IR  = 
1.94  dB;  (b)  S3IR  =  1.87  dB;  (c)  S3IR  =  1.67  dB. 
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dpuetc,x)  lipue  b  tcdpue 


a2  3  =  0.8  ,  S3I2R=  1.94  dB  =  0.4,  S^R-  7.96  dB 


a2,3  =  0-8’S3'2R=1-94dB 


a  =  0.6,  SIR  =  4.44  dB 

1,3  3  1 


a2 3  =  0.8  ,  S3I2R  =  1 .94  dB  a,  3  =  0.8,  S^R  =  1 .94  dB 


Figure  22.  Plots  of  the  actual  gain  parameters  a2  3 ,  al3  and  the  estimators  /?2  3 ,  /?, ,  vs. 
different  S3NR  values  for  multiple  simulations  of  y3  while  a23  =  0.8:  (a)  S3IR  = 
0.97  dB;  (b)  S3IR  =  0  dB;  (c)  S3IR  =  -1.07  dB. 
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Figure  23. 


(b| 

Variance  of  estimators  J323,J313  vs.  their  respective  SNR  with  a2  3  =  0.2,  S3I2R 
=  13.98  dB:  (a)  a13  =  0.1,  S3I1R  =  20  dB;  (b)  al3  =  0.4,  S3I1R  =  7.96  dB. 


<b) 

Figure  24.  Variance  of  estimators  J323,J313  vs.  their  respective  SNR  with  a2  3  =  0.5, 
S3I2R  =  6.02  dB:  (a)  al  3  =  0.1,  S3I1R  =  20  dB;  (b)aI  3  =  0.4,  S3I1R  =  7.96  dB. 
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A  summary  of  SER  performances  for  multiple  configurations,  each  with  different 
combinations  of  gain  parameter  values  is  shown  in  Figure  25.  The  common  quantity  in 
the  plots  is  the  S3IR  experienced  by  each  of  these  configurations.  The  best  SER  curves 
correspond  to  high  S3IR  values,  which  points  out  the  fact  that  in  order  to  maintain  a 
system  with  accurate  estimations,  a  reasonable  energy  difference  is  needed  between  the 
main  signal  and  the  interfering  ones.  Specifically,  configurations  with  S3IRS  values  of 
6.99  dB  or  higher  have  SER  curves  which  follow  the  theoretical  curve  quite  well. 


Figure  25.  Performance  curves  of  SER  vs.  S3NR  for  post  s3  decisions  at  multiple 
configurations:  (a)  a23  =  0.01;  (b)  a23  =  0.2;  (c)  a2  3=0.5;  (d)  <x,3=0.8. 
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IV.  ADDITION  OF  A  FOURTH  RECEIVER 


A.  PARAMETER  ESTIMATION  FOR  A  FOUR-RECEIVER  SYSTEM 

Here  we  consider  that  a  fourth  receiver  is  added  to  the  multi-platform  system  and 
takes  a  position  somewhere  in  Sector  D  from  Figure  26  (shown  again  for  convenience.) 


Figure  26.  Configuration  of  multiple  BSs  with  overlapping  antenna  coverage  which  has 
different  signal  combinations  received  in  each  sector.  From  [3]. 


If  the  received  signal  (2.4)  is  taken,  and  the  LSE  method  is  applied,  we  get 


^(«)=||y4-(S4+«3,4S3+«2,4S2+«l,4Sl,||  (4  ^ 

(y4  ^4  ^3,4^3  ^2,4^2  ^1,4^1 )  (y 4  ^4  ^3,4^3  ^2,4^2  ^I.Al  )' 


Despite  a  lengthy  derivation,  it  is  shown  that  the  estimators  for  the  gain  parameters  are 


^3,4 


Re{y/s3 }  -  Re{/?4-3 }  -  a24  Re{  a  2 }  -  a,4  Re  (a,,  } 


'4,3  J  **2,4 

2  E. 


oc2A 


Re{y/s2}-Re{p42}-«34Re{A,2}-«i,4  Re(P2,i} 

2  E. 


(4.2) 


ai,4 


Re{y4Hs1 }  - Re{  A 1 }  - a3  4  Re{p31  }-a24  Re{p2l } 


2  Es 
S1 
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Unlike  in  Equation  (3.8),  which  is  solved  via  simple  substitution.  Equation  (4.2)  does  not 
lend  itself  to  simple  derivation.  Conveniently,  the  three  equations  and  three  unknowns 
allow  us  to  derive  the  estimates  via  linear  algebraic  method.  By  reordering  (4.2)  to  move 
the  gain  parameters  to  one  side,  we  get 

Re{y/s3 }  -  Re{/?4 3 }  =  a3A(2ES} )  +  a24  Re{p3  2 }  +  al4  Re{p3 4 } , 

Re{y4"s2 }  -  Re{/?4  2 }  =  a3A  Re{ /?3  2 }  +  a24(2ESi )  +  al  4  Re{/?2 4 } ,  (4.3) 

Re{y4"si }  - Re{  A.i )  =  «3,4  Re{ A.i  1  +  aiA  Re( A.i  1  +  aiA(2ESl  )• 

For  simplicity,  (4.3)  is  arranged  to 

~A~|  |"X,  Fj  Z,  1  a3A 

B  =  F,  Z2  a24  ,  (4.4) 

C  X3  Y3  Z3  a14 

where  we  define  A  =  Re{y4Hs3}-Re{/?43},  5  =  Re{y4//s2}-Re{p42}, 

C  =  Re{y/s1}-Re{/?41},  Xx=2 ESs,  Yt=Re{p32},  Z,  =  Re{Al},  X2=Re{A2}, 

Y2=2Es ,  Z2  =  Re{/?21} ,  X3  =  Re{yo31}  ,  F3  =  Re{p21} ,  and  Z3=2E^  .  To  successfully 

solve  for  the  three  unknowns,  the  inverse  of  the  X,  Y,  and  Z  matrix  of  constants  is 
multiplied  with  the  A,  B,  C  vector  of  constants  such  that 

a3A  rx,  y{  zx  rpr 

a2,4  =  *2  K  Z2  B  .  (4.5) 

a{  4  X3  F3  Z3  _  C 

Similar  to  Chapters  II  and  III,  there  is  a  need  for  Rx4  to  make  a  decision  on  the  received 
signal  y4  and  create  a  pre  calculation  for  s4 ,  or 

s4  =  dec(y4),  (4.6) 

which  when  substituted  into  (4.5)  yields  the  estimates 
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>3,4 1  IX  Y,  Z|  IXa- 

A,4  =  z2  r2  Z2  B  ,  (4.7) 

A,4  _  X3Y3Z3  c 

where  A  =  Re{y/s3 }  -  Re{p43 } ,  B  =  Re{y4% }  -  Re{p4  2 } ,  C  =  Re{y/s, }  -  Re{p4 , } , 

Xx  =  2ES 3,  Y{  =Re{/332},  Z,=Re{p31},  -Y2  =  Re{/?32}  ,  Y2=2ES^,  Z2  =  Re{/)21}, 

X3  =  Re{/?3  j } ,  F3  =  Re{/321} ,  and  Z3=2ES  .  The  solutions  from  (4.7)  are  now  used  in  the 
computation  of  the  post  decisions,  or 

s4  =  dec(y4  -j334 s3  - j324s2  -/?li4Sj).  (4.8) 

B.  RESULTS  FOR  A  FOUR-RECEIVER  SYSTEM 

Following  the  same  simulation  setup  from  Chapter  III,  we  assume  that  Rx4  is 
provided  the  reference  signals  s3 ,  s2 ,  s,  that  are  mostly  interference-free. 


1.  Fix  S4NR,  Vary  a3  4,  a2,4  and  <*,4 

The  same  approach  from  Chapter  III  is  used  to  calculate  the  total  SIR,  such  that 


S4IR  = 


E.  +E.  +E. 


1  ^3,4^3  "t"  ^2,4^2  ^1,4^1 


With  the  energy  for  each  signal  set  to  the  same  value,  (4.9)  becomes 


S4IR 


(4.10) 


Again,  by  first  setting  specific  noise  levels  in  the  simulation,  we  fix  S4NR  values 
and  compare  the  accuracy  of  the  estimates  against  theoretical  values.  From  Figure  27,  we 
see  /?34  calculations  versus  S4I3R  values,  and  that  at  high  S4NR  values,  the  estimations 

only  show  good  accuracy  when  the  energy  of  s4  is  much  greater  than  the  interference  in 
the  rest  of  the  system  (i.e.,  low  a  values).  The  performance  of  the  estimators  show 
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increased  degradation  when  more  noise  is  added  to  the  system  and  the  S4NR  values  are 
reduced.  This  is  shown  in  Figure  28.  Taking  a  look  at  the  direct  comparision  between  the 
calculations  of  fi34  and  the  actual  values  of  a34,  we  see  similar  results  wherein  high 

S4NR  only  provides  accurate  estimates  around  gain  parameters  values  of  0.5  as  seen  in 
Figure  29.  Likewise,  increased  noise  in  the  system  makes  accurate  gain  estimation  very 
difficult  as  calculations  made  when  there  is  low  S4NR  diverge  almost  immediately  from 
the  theortical  line  in  Figure  30. 

If  a2  4  is  increased  to  0.4,  the  overall  interference  in  the  system  increases  and  we 

again  see  degradation  in  the  accuracies  of  each  estimate.  As  shown  in  Figures  31-34,  the 
three  interfering  signals  demonstrate  their  ability  to  drastically  reduce  the  performance  of 
the  estimators  by  compounding  the  effects  of  the  amplitude  gain  parameters.  Even  at  high 
S4NR  levels,  the  estimates  barely  approach  the  a34  line  unless  the  actual  gain  is  set  very 

low  to  begin  with  in  the  simulation. 

It  is  important  to  emphasize  the  fact  that  because  of  the  relationship  examined  in 
(4.10),  fixed  SNR  simulations  do  not  require  separate  runs  for  each  of  the  individual  gain 
parameters.  Due  to  the  overall  SIR  caused  by  the  additive  effects  of  the  interference 
gains,  holding  a24  and  a14  at  certain  values  while  cycling  through  several  a34  values 

would  provide  similar  (if  not  exact)  results  as  cycling  through  a  different  gain  parameter 
and  holding  the  others  constant.  These  connections  are  further  examined  in  the  Appendix. 
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(a) 


S4l3R|dB| 

|c| 


S4l3R(dB) 

(b| 


S4l3R(dB) 

(d) 


Figure  27.  Plots  of  gain  estimates  j334  vs.  S4I3R  values  with  fixed  high  S4NR  levels  while 
a24  =  0.2:  (a)  alA=  0.01;  (b)a14=  0.  1;  (c)a14=  0.2;  (d)«14=0.4. 


S4I3R  fdB) 
(b| 


S4I3R  (dB) 
(d) 


Figure  28.  Plots  of  gain  estimates  J334  vs.  S4I3R  values  with  fixed  low  S4NR  levels  while 
a2 4  =  0.2:  (a)  alA=  0.01;  (b)  al4=  0.  1;  (c)  al4=  0.2;  (d)  al4=  0.4. 
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“3,4 

(a) 


“3,4 

(c) 


“3,4 

(b) 


“3,4 

(b) 


Figure  29.  Plots  of  gain  estimates  [3:,  4  vs.  actual  a  ,  4  values  with  fixed  high  S4NR  levels 
while  a2A  =  0.2:  (a)  al4  =  0.01;  (b)  al4  =  0.1;  (c)  al4  =  0.2;  (d)  al4  =  0.4. 


“3.4 

(a) 


“3.4 

(b| 


“3,4  “3,4 

(C)  (d) 


Figure  30.  Plots  of  gain  estimates  4  vs.  actual  a3  4  values  with  fixed  low  S4NR  levels 
while  a2  4  =  0.2:  (a)  a14  =  0.01;  (b)  al4  =  0.1;  (c)  al4  =  0.2;  (d)  al  4  =  0.4. 
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(a) 


S4l3R(dB) 

<b) 


S4l3R(dB) 

(d) 


Figure  31.  Plots  of  gain  estimates  j334  vs.  S4I3R  values  with  fixed  high  S4NR  levels  while 
a2A  =  0.4:  (a)  alA  =  0.01;  (b)  a14  =  0.  1;  (c)  alA=  0.2;  (d)  alA=  0.4. 


S4I3R  (dB) 
(a) 


Figure  32.  Plots  of  gain  estimates  j33  4  vs.  S4I3R  values  with  fixed  low  S4NR  levels  while 
a24  =  0.4:  (a)  al4=  0.01;  (b)  al4=  0.  1;  (c)  alA=  0.2;  (d)  al4=  0.4. 
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Figure  33. 


“3,4  “3.4 

(c|  (d) 

Plots  of  gain  estimates  ft3  4  vs.  actual  a3  4  values  with  fixed  high  S4NR  levels 
while  a24=  0.4:  (a)  alA  =  0.01;  (b)  al4  =0.1;  (c)  al4  =  0.2;  (d)  al4  =  0.4. 


“3,4 

(a) 


“3.4 

(b) 


“3.4  “3.4 

(c)  (d) 

Figure  34.  Plots  of  gain  estimates  ft3  4  vs.  actual  a3  4  values  with  fixed  low  S4NR  levels 
while  a2A  =  0.4:  (a)  a14  =  0.01;  (b)  al4  =  0.1;  (c)  al  4  =  0.2;  (d)  al  4  =  0.4. 
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2. 


Fix  «3  4 ,  0.2,4  and  Oj  4 ,  Vary  S4NR 


Observing  the  SER  of  the  post  calculations  for  each  of  the  combinations  of  gain 
parameter  values,  we  better  grasp  the  severity  of  the  compounding  effects  of  interference 
on  the  SOL  Only  at  high  S4IR  levels  (i.e.,  greater  than  6.0  dB)  do  the  generated  s4  signals 

provide  results  that  are  considered  acceptable  for  use.  For  the  four-receiver  system,  it  is 
unnecessary  to  examine  the  pre  calculation  results  due  to  expected  poor  performance  in 
such  an  interference-filled  environment.  The  simulations  that  are  conducted  increase  one 
or  two  gain  parameters  at  a  time,  effectively  decreasing  the  accuracy  of  each  of  the 
estimators  of  Equation  (4.7)  and  therefore  the  performance  of  SER  is  also  decreased. 
From  Figure  35,  it  is  seen  that  at  relatively  low  a  values,  the  SER  for  the  post 
calculations  follow  the  theoretical  line  for  a  QPSK  system  with  no  interference.  Once  the 
S4IR  drops  below  6.0  dB,  the  performance  starts  to  degrade  considerably.  Increasing  the 
total  interference  in  the  system,  we  see  that  the  SER  continues  to  have  a  steady  decrease 
and  shifts  away  from  an  interference-free  system,  as  seen  in  Figures  36  and  37. 


Figure  35.  Performance  curves  of  SER  vs.  S4NR  for  post  decision  of  s4  at  different 
interference  combinations  when  a34  =  0.1:  (a)  a24  =  0.2;  (b)  a24  =  0.5. 
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SER  SER 


S4NR  (dB) 

(a) 


OL  .3 

LU  10 
OT 


S4NR  (dB) 
(b) 


Figure  36.  Performance  curves  of  SER  vs.  S4NR  for  post  decision  of  s4  at  different 
interference  combinations  when  a34  =  0.3:  (a)  a24  =  0.2;  (b)  a2  4  =  0.5. 


—  —  -  No  Interference 

—•—a.,  ,  =  0.1,  S,IR  =  2.08  dB 

1.4  4 

—*—<2.,  ,  =  0.4,  SIR  =  1.14  dB 
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Figure  37.  Performance  curves  of  SER  vs.  S4NR  for  post  decision  of  s4  at  different 
interference  combinations  when  a34  =  0.6:  (a)  a2  4  =  0.2;  (b)  a24  =  0.5. 
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Next  we  observe  the  accuracy  of  the  estimations  against  the  actual  gain 
parameters  that  is  used  in  each  simulation.  The  same  collection  of  gain  values  from  the 
SER  performance  evaluations  are  used  to  show  situations  with  both  low  and  high 
interference.  Each  row  from  Figures  38-43  represents  a  different  configuration  of  the 
simulation.  For  example,  from  Figure  38  we  see  two  instances  of  the  modeled  system 
with  the  top  row  showing  gain  values  and  estimates  for  a34  =0.1,  a24  =  0.2  and  al4  = 

0.1  while  the  bottom  row  shows  a  different  run  with  a34  =0.1,  a,  4  =0.2  and  ax  4  =  0.4. 

The  performance  of  the  estimations  generally  show  increased  accuracy  when 
higher  S4NR  or  S4IR  is  used.  As  expected,  once  the  S4IR  drops  below  6.0  dB,  the  quality 
of  the  estimations  begin  to  deteriorate.  Some  configurations  even  experience  a  reduction 
of  accuracy  with  the  increase  of  S4NR  as  seen  in  Figure  43.  Again,  this  is  likely  due  to 
the  inability  of  the  estimators  to  handle  scenarios  where  the  energy  of  the  SOI  does  not 
strongly  overpower  the  other  interference  and  noise  energies.  It  should  be  remembered  as 
well  that  the  actual  gain  parameter  estimations  fl  are  calculated  with  use  of  (4.6),  which 
is  heavily  interfered  by  three  other  signals  other  than  s4 .  In  the  next  section,  we  discuss 
another  technique  to  take  advantage  of  the  cleaner  s4  calculation  in  order  to  recalculate  a 
better  gain  estimate. 
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Figure  38.  Plots  of  the  actual  gain  parameters  a34,  a3  4  ,  al4  and  the  estimators  ft,  4 ,  fl2  4  ,  /?,  4  vs.  different  S4NR  values  for  multiple 

simulations  of  y4 :  (a)  S4IR  =  12.2  dB;  (b)  S4IR  =  6.78  dB. 
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Figure  39.  Plots  of  the  actual  gain  parameters  a34,  a2  4  ,  a14  and  the  estimators  ft,  4 ,  //,  4  ,  J3l4  vs.  different  S4NR  values  for  multiple 


simulations  of  y4 :  (a)  S4IR  =  5.69  dB;  (b)  S4IR  =  3.77  dB. 
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Figure  40.  Plots  of  the  actual  gain  parameters  a34,  a2  4  ,  al4  and  the  estimators  fi,  4 ,  [i,  4  ,  /?,  4  vs.  different  S4NR  values  for  multiple 


simulations  of  y4 :  (a)  S4IR  =  8.54  dB;  (b)  S4IR  =  5.38  dB. 
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Figure  41 .  Plots  of  the  actual  gain  parameters  a3  4 ,  a2  4  ,  a14  and  the  estimators  ft,  4 ,  fi,  4  ,  /?,  4  vs.  different  S4NR  values  for  multiple 


simulations  of  y4 :  (a)  S4IR  =  4.56  dB;  (b)  S4IR  =  3.01  dB. 
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Figure  42.  Plots  of  the  actual  gain  parameters  a34,  a2  4  ,  al4  and  the  estimators  fi,  4 ,  /?2  4  ,  (3l4  vs.  different  S4NR  values  for  multiple 


simulations  of  y4 :  (a)  S4IR  =  3.87  dB;  (b)  S4IR  =  2.52  dB. 
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Figure  43.  Plots  of  the  actual  gain  parameters  a34,  a24,  al4  and  the  estimators  ft,  4 ,  /32  4  ,  /i,  4  vs.  different  S4NR  values  for  multiple 

simulations  of  y4 :  (a)  S4IR  =  2.08  dB;  (b)  S4IR  =  1.14  dB. 
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C.  RESULTS  OF  ITERATIVE  ESTIMATION  CALCULATION 

With  the  combination  of  the  multiple  interferences  negatively  affecting  the 
performance  of  the  estimators,  another  technique  is  employed  to  improve  the  accuracy  of 
the  /?  calculations.  From  Equation  (4.7),  we  see  that  the  estimators  make  use  of  the  pre 
calculation  s4  in  order  to  make  a  calculation.  By  inspection,  it  is  clear  this  is  not  optimal 
as  Rx4  is  required  to  make  a  decision  based  simply  on  the  value  of  y4  which  is  heavily 
interfered  by  the  other  signals  in  the  environment.  Therefore,  by  reusing  the  cleaner  post 
decision  s4  with  the  estimate  equations,  we  are  able  to  produce  a  more  accurate  estimator 

through  a  second  round  of  calculation.  Substituting  s4  with  Equation  (4.8)  and 
incorporating  it  back  into  (4.7),  we  get 

PlA 
PlA 

PlA 

where  the  updated  A,  B  and  C  terms  are  now  4  =  Re  { v }  -  Re  { p43 } , 

5  =  Re{y4ws2}-Re{/34  2}  and  C  =  Re{y4ws1}-Re{/?41}.  Increasing  the  number  of 

recalculations,  or  rounds,  also  continues  to  improve  performance  of  the  estimators,  in 
turn  creating  even  cleaner  estimates  of  s4 .  The  Monte  Carlo  simulations  are  performed 

again,  this  time  iterating  the  decision  process,  providing  updated  estimate  and  re¬ 
demodulated  post  signals.  As  can  be  seen  in  Figures  44  and  45,  two  Pn  levels  are  set  to 
illustrate  one  system  with  high  SNR  and  another  with  low  SNR.  The  first  and  third 
iterated  rounds  are  illustrated  and  it  is  clear  that  the  use  of  s4  in  (4.11)  vastly  improves 
the  quality  of  the  estimates  in  Round  3  (Rd3). 

Even  with  a2  4  =  0.7,  the  performance  of  the  estimators  at  the  end  of  Rd3  shows 

the  considerable  improvement  in  estimation  accuracy  as  seen  in  Figures  46  and  47. 
Without  this  iterative  recalculation,  there  is  little  chance  of  successfully  extracting  SOIs 
from  the  received  signals. 
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Figure  44.  Iterative  gain  estimates  of  )3,  4  vs.  a3  4  with  fixed  S4NR  levels  while 
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Figure  45.  Iterative  gain  estimates  of  j334  vs.  S4I3R  with  fixed  S4NR  while  a24  =  0.4: 

(a)  a, ,  =  0.2;  (b)  a, ,  =  0.4. 
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Figure  46.  Iterative  gain  estimates  of  /7,  4  vs.  cr3  4  with  fixed  S4NR  levels  while 
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Figure  47.  Iterative  gain  estimates  of  /?3  4  vs.  S4I3R  with  fixed  S4NR  while  a1 4  =  0.7: 

(a)  alA  =  0.2;  (b)  aM  =  0.4. 
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To  really  see  the  positive  effect  of  this  iterative  gain  estimation,  we  again  view 
the  estimates  as  they  are  recalculated  alongside  the  actual  gain  parameter  used  in  the 
specific  configuration.  As  seen  in  Figures  48-51,  each  row  demonstrates  the  three 
amplitude  gain  parameters  associated  with  the  interfering  signals  of  y4  for  one  run  of  the 

simulation.  Each  plot  also  shows  the  results  of  the  original  gain  estimate  calculation  as 
well  as  two  additional  rounds  (i.e.,  Rd2  and  Rd3)  of  re-estimation.  For  the  configurations 
which  previously  suffered  immense  performance  degradation  due  to  high  amount  of 
interference  in  the  received  signal,  this  new  technique  greatly  enhances  estimation 
accuracy.  Even  when  the  S4IR  is  very  low  (e.g.,  Figure  51b  with  0.46  dB),  the  third  round 
of  iterative  estimate  calculation  provides  nearly  perfect  estimates  at  high  S4NR. 

In  addition  to  the  improved  accuracy  of  the  gain  estimates,  another  result  from 
this  iterative  technique  is  that  each  new  round  also  yields  a  more  accurate  s4  calculation. 

The  performance  curves  for  this  IC  approach  show  significantly  improved  SER  for  low  to 
medium  values  of  S4NR  levels,  as  seen  in  Figures  52  and  53.  The  resulting  SERs  from 
the  second  and  third  rounds  now  follow  the  theoretical  curve  from  low  to  medium  S4NR. 
However,  it  is  important  to  notice  that  at  high  S4NR  levels,  there  appears  to  be  a 
performance  “floor”,  where  no  amount  of  reduction  of  noise  reduces  SER  below  about 
10°.  We  believe  that  is  likely  due  to  the  fact  that  sp  s2,  s3  are  themselves  estimates 
which  are  slightly  corrupted  with  specific  noise  power  in  the  simulation.  Perhaps  if  these 
reference  signals  were  set  up  with  higher  SNR,  then  the  performance  floor  could  be 
lowered. 
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Figure  48.  Plots  of  the  actual  gain  parameters  a34,  a24,  a14  and  the  iterated  estimations  of  /?3  4 ,  /324,  /3l4  vs.  different  S4NR 

values  for  multiple  simulations  of  y4 :  (a)  S4IR  =  6.78  dB;  (b)  S4IR  =  4.44  dB. 
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Figure  49.  Plots  of  the  actual  gain  parameters  a34,  a24,  al4  and  the  iterated  estimations  of  ft3  4 ,  ft24,  ftl4  vs.  different  S4NR  values 

for  multiple  simulations  of  y4 :  (a)  S4IR  =  3.77  dB;  (b)  S4IR  =  2.44  dB. 
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Figure  50.  Plots  of  the  actual  gain  parameters  a34,  a24,  al4  and  the  iterated  estimations  of  4 ,  /?,  4  ,  4  vs.  different  S4NR  values 

for  multiple  simulations  of  y4 :  (a)  S4IR  =  2.68  dB;  (b)  S4IR  =  1.61  dB. 
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Figure  51.  Plots  of  the  actual  gain  parameters  a34,  a24,  al4  and  the  iterated  estimations  of  /?3  4 ,  f32  4,  J3l4  vs.  different  S4NR  values 

for  multiple  simulations  of  y4 :  (a)  S4IR  =  1.25  dB;  (b)  S4IR  =  0.458  dB. 
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SER 


Figure  52.  Performance  curves  of  SER  vs.  S4NR  for  iterated  s4  calculation  at  different 
interference  combinations  when  a34  =  0.4:  (a)  «24  =  0.2;  (b)  a24  =  0.5. 


Figure  53.  Performance  curves  of  SER  vs.  S4NR  for  iterated  s4  calculations  at  different 
interference  combinations  when  a34  =  0.7:  (a)  a24  =  0.2;  (b)  a24  =  0.5. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  success  of  RSIC  relies  on  the  initial  reception  of  an  interference-free 
reference  signal  and  the  accuracy  of  the  amplitude  gain  estimation  of  the  interfering 
signal.  This  estimation  is  crucial  in  correctly  subtracting  the  interference  a  system 
experiences  from  multiple  transmitters.  It  was  observed  that  factors  such  as  noise  and 
compounding  interference  greatly  impact  the  accuracy  of  the  estimators  used  in  the 
Monte  Carlo  simulations. 

A.  CONCLUSIONS 

For  this  thesis,  it  was  assumed  that  clean  reference  signals  were  available  for 
subtracting  interference  from  subsequent  receivers.  Obviously,  the  placement  of  said 
receivers  in  favorable  locations  is  vital  to  ensuring  cancellation  of  wanted  interference. 
For  the  best  results,  it  is  advised  to  have  the  receiver  placed  so  that  the  energy  of  the 
desired  signal  is  much  greater  than  the  energy  of  interfering  signals  in  the  area. 

The  LSE  estimation  method  was  utilized  in  this  thesis  and  while  capable  of 
finding  accurate  estimates  for  systems  with  low  interference,  its  limitation  was  observed 
when  attempting  to  estimate  values  of  oc  that  approached  1.0.  At  these  levels,  the 
amount  of  interference  in  the  received  signals  affected  the  symbols  such  that  calculations 
for  the  estimators  could  not  produce  a  satisfactory  (3  value.  Introduction  of  an  iterative 
gain  estimation  technique  vastly  improved  the  performance  of  estimates  in  a  four- 
receiver  system.  The  technique  may  also  work  for  whatever  number  of  receivers. 

Cleary,  the  noise  and  interference  present  in  the  systems  play  a  major  part  in 
determining  the  amplitude  gain  parameters,  and  the  higher  the  SNR  values  are  for  a 
signal,  the  more  accurate  an  estimation  is  expected  to  become.  In  turn,  this  provides  for  a 
more  accurate  post  demodulation  the  desired  signal. 

B.  RECOMMENDATIONS 

While  QPSK  is  still  a  relevant  modulation,  more  powerful  and  complex 
communication  methods  exist  that  would  most  likely  perform  differently  with  the  IC 
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technique  described  in  this  thesis.  Follow  on  work  that  focuses  on  SIC  for  multi-carrier 
and  MIMO  systems  would  be  most  beneficial. 

Other  estimation  methods  exist  that  could  possibly  provide  better  estimator 
performance  for  this  multi-platform  system.  Future  work  could  observe  the 
implementation  of  different  methods  to  better  estimate  a  ,  perhaps  at  levels  higher  than 
1.0,  or  when  the  gain  is  complex. 
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APPENDIX 


Due  to  the  relationship  of  the  terms  in  Equation  (3.16)  for  a  three -receiver  system, 
certain  combinations  of  gain  parameters  yield  similar  if  not  exact  performance  SER 
outcomes  as  seen  in  Figure  54. 


(c)  (a) 

Figure  54.  Performance  results  observed  from  Figure  11c  with  a23  =  0.01  and  a13  =  0.2 
compared  against  Figure  12a  where  a2  3  =  0.2  and  at  3  =0.1.  The  S3IR  for  both 

simulations  is  14  dB. 

The  same  outcome  is  illustrated  in  terms  of  the  performance  of  the  estimators, 
where  switching  the  values  between  a2  3  and  at  3  would  result  in  related  accuracies  of  the 

gain  estimations,  as  is  seen  in  Figure  55.  This  association  mitigates  the  requirement  of 
multiple  simulations  with  each  gain  parameter  individually  being  held  constant. 
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dpUBttn  t6dpUB 


a2  3  =  0.01 ,  S3I2R  =  40  dB  a,  3  =  0.8,  S^R  =  1 .94  dB 


a2  3  =  0.8  ,  S3I2R  =  1 .94  dB  a13  =  0.01 ,  S^R  =  40  dB 


Figure  55.  Performance  results  observed  from  Figure  16c  with  a2  3  =  0.01  and  al3  =  0.8 
compared  against  Figure  21a  where  a2  3  =  0.8  and  3  =  0.01 .  The  S3IR  for  both 

simulations  is  1.94  dB. 


In  a  four-receiver  system,  this  same  relationship  is  observed  in  Equation  (4.10), 
and  the  performances  of  the  SER  and  estimators  are  seen  in  Figures  56  and  57. 


72 


SER 


Figure  56.  Performance  results  observed  from  Figure  36a  with  a34  =  0.3,  a24  =  0.2  and 
compared  against  Figure  56b  (a  separate  configuration  simulation)  where  a34  = 
0.2  and  a24  =  0.3.  The  S4IRS  for  both  simulations  are  8.54  and  5.38  dB  for  a  14  = 

0.1  and  0.4,  respectively. 
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ct  .  =  0.3,  SIR  =  10.5  dB  a  ,  =  0.2  ,  SIR  =  14  dB 
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a,  =0.1  ,  S.I.R  =  20  dB 

1,4  4  1 


a  ,  =  0.2,  SIR  =  14  dB  a  =  0.3  ,  SIR  =  10.5  dB 

3,4  7  4  3  2,4  7  4  2 


a  =  0.1  ,S  ,1  R  =  20  dB 

1,4  ’41 


(b) 

Figure  57.  Performance  results  observed  from  Figure  40a  with  a34  =  0.3,  a24  =  0.2  and  al4  =0.1  compared  against  Figure  57b  (a 


separate  configuration  simulation)  where  a3  4  =  0.2,  a24  =  0.3  and  al4  =  0.1.  The  S4IR  for  both  simulations  is  8.54  dB. 
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